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Strangeness enhancement in heavy-ion collisions is stud- 
ied at the parton level by examining the partition of the new 
sea quarks generated by gluon conversion into the strange and 
non-strange sectors. The CTEQ parton distribution functions 
are used as a baseline for the quiescent sea before gluon con- 
version. By quark counting simple constraints are placed on 
the hadron yields in different channels. The experimental val- 
ues of particle ratios are fitted to determine the strangeness 
enhancement factor. A quantitative measure of Pauli blocking 
is determined. Energy dependence between SPS and RH1C 
energies is well described. No thermal equilibrium or statisti- 
cal model is assumed. 

PACS numbers: 25.75.Dw, 24.85. +p 

I. INTRODUCTION 

The production of strange particles in heavy-ion col- 
lisions has been a subject of intense study in the past 
twenty years, ever since the proposal that it may reveal 
a signal of quark-gluon plasma formation (l],^]. Various 
approaches to the problem have been adopted, ranging 
from statistical thermal model j^,^] to simple quark co- 
alescence model HH to dual parton model JtJ. Despite 
differences in diverse viewpoints, the major theme is to 
explain the phenomenon of strangeness enhancement in 
nuclear collisions Although the experimental defini- 
tion of strangeness enhancement is the increase of the 
strangeness content of the produced hadrons with in- 
creasing number of participants from pp to A A collisions, 
a more appropriate theoretical description of strangeness 
enhancement is in terms of the increase of strange quarks 
before hadronization. In this paper we present a quanti- 
tative treatment of the enhancement factor in the frame- 
work of the parton model, and obtain a numerical mea- 
sure of Pauli blocking. 

The point stressed in the original explanation for 
strangeness enhancement [^J|2j is that when the quark 
degrees of freedom are liberated, it is easier to create 
strange quark pairs than strange hadrons because the 
ss threshold is lower. Deconfincmcnt then leads natu- 
rally to the possibility of plasma formation. In our view 
the quarks have always been the basis for understand- 
ing hadron production even in pp collisions for ^/s > 10 
GeV. The recombination model has been able to repro- 
duce the \ow-px inclusive distributions in the fragmenta- 



tion region by treating the hadronization processes at the 
parton level . Thus the relevance of the quark de- 
grees of freedom in heavy-ion collisions at SPS is nothing 
new. In collisions at such energies the nucleons are bro- 
ken up and thus deconfined, but it does not mean that 
there is thermalized quark-gluon plasma, which no one 
would associate with pp collisions. 

Once one descends to the parton level, the notion of 
strangeness enhancement (SE) can take on a quantita- 
tive description in terms of the strange quark population. 
The baseline for the unenhanced s quark distribution in 
the quiescent sea should be pinned down by the parton 
distribution functions of the nucleon studied exhaustively 
by several groups |T^-|l^]. We shall use the distribution 
functions of the CTEQ global analysis that fits some 
1300 data points obtained for many reactions in 16 ex- 
periments. Their extrapolations to Q 2 — 1 (GeV/c) 2 
are presented in the form of graphs available on the web 
]l4]| . Since gluons do not hadronize directly, there be- 
ing no glueballs found, they are converted to quark pairs 
which subsequently hadronize by recombination. How 
much the conversion goes into the strange sector gives us 
a measure of SE. 

Gluon conversion is not a new process that we must 
consider for heavy-ion collisions. Even in hadronic col- 
lisions gluons must convert in order to hadronize. Such 
conversion has been included in the study of inclusive dis- 
tributions of hadrons produced in the fragmentation re- 
gion in the framework of the valon-recombination model 
pOI , and more recently using the CTEQ parton distri- 
bution functions JlJ] to reproduce various hadronic spec- 
tra |l5j ]. Our attention in this paper is shifted from the 
fragmentation region in hadronic collisions, where the x 
dependence is an issue, to the central region in nuclear 
collisions, where the relative yields of particles produced 
are the focus. 

Clearly, there is no way to study SE from first prin- 
ciples. Our investigation here is phenomenological. Our 
goal is very modest. It is not possible to compute parti- 
cle ratios from the parton model alone. We shall use a 
large set of particle ratios as experimental inputs to guide 
us in the determination of the SE factor. The effect of 
Pauli blocking in the non-strange sector is included in the 
those inputs, and are not amenable to first-principle cal- 
culations. Our theoretical input is essentially the count- 
ing of quarks and antiquarks in their partition into the 
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various hadronic channels. In that sense the physics in- 
volved is basically the same as in the coalescence model 
which is a simplified version of the recombination 
model S. The emphases in Refs. (H)^] are in the multi- 
plicative aspect of the probabilities of having quarks and 
antiquarks in the same region of phase space in their for- 
mation of hadrons. That results in an undesirable feature 
of s and s imbalance in the linear version j^] , which is not 
satisfactorily resolved in the nonlinear version q] by the 
introduction of unknown factors. Our emphasis here is 
on the partition of the q,q, s,s quarks into the hadronic 
channels and on the enhancement of their populations 
from gluon conversion. We shall not investigate the im- 
plications of the hadron probabilities being products of 
the quark probabilities. 



II. QUARK COUNTING 

We begin by drawing the boundary of our concern here. 
Since the yields on multi-strange hyperons are low com- 
pared to K and A, we shall in first approximation ignore 
the production of S and f2, and aim at results with ac- 
curacies not better than 90%. With such simplification 
we can better exhibit the spirit of our approach to the 
problem and make more transparent the issues involved 
in SE. Improvements that include the S and Q particles 
can be considered later. We shall also consider isosym- 
metric dense medium at mid-rapidity so that we need not 
distinguish u and d quarks. Proton and neutron will be 
equal in number, as do 7r + and ir~ . The strange quarks 
s and s are produced in equal numbers, but K + and K~ 
will not be produced in equal numbers because of asso- 
ciated production. 

Let us use the following notation to denote the num- 
bers of hadrons and non-strange quarks, e.g., N is the 
number of nucleons, and q is the number of light quarks. 



N = p + n, 

II = 7T + + 7T~ + 7T 

Y = A + E° + £+ 
K = K+ + K°, 
q = u + d. 



N 



P- 



K = K~ 



(1) 
(2) 
(3) 
(4) 
(5) 



Then there are linear relations among these numbers 
based on counting the number of valence quarks in the 
various hadrons 



3N + n + 2Y + K = q, 
3N + IT + 2Y + K = q, 

Y + K = s 

Y + K = s. 



(6) 
(7) 
(8) 
(9) 



The right-hand sides of the above equations all refer to 
the numbers of quarks after enhancement from gluon con- 
version. Let k be the fraction of s quarks that recombine 



with non-strange antiquarks to form anti-kaons, and sim- 
ilarly R be the fraction of s to form kaons. That is, we 
define 



K = k s, K = k s. 

Then on account of Eqs. (8) and (9), we have 

Y = (1 - k) s, Y=(l-R)a. 

Define the hadronic ratios 

r = K/K, R = Y/Y. 

It then follows that 
1 - R 



K 



R 



k = r k, 



(10) 

(11) 
(12) 

(13) 



where s = s has been used. For experimental values of 
r and R, we assume that K/K sw K + / K~ and Y/Y 
A/A. For Pb-Pb collisions at SPS the values are ||-[l|] 

r = 1.8, R = 0.13, (14) 

so we obtain 

K = 0.52, R = 0.94. (15) 

With these values of k and R we can proceed to consider 
the non-strange sector. Define 



p = N/N 

so that we can obtain from Eqs. (6) and (7) 
1 



N 

n 



3(l-p) 

1 t 

{-pq + q 



[q — q + 3 s (k — k)] , 



1 



P 



-s [(1 + 2 p) K - (2 + p) R + 2(1 - p)}} . 



(16) 



(17) 



(18) 



The particle ratios involving abundant strange and non- 
strange hadrons are 



1 



1-P 
1 



N _ 
K 
II 

K ~ (1 - p)R 



3sR 



1 

- - 1 

r 



(19) 



(1-P) 



1 <b> 



s s 

- (l + 2p)/s+(2 + p)/s-2(l-p)], (20) 

Where q v denotes the number of valence quarks, i.e., 
q v = q — q. The RHS can be determined from parton 
distributions, assuming that the parameters k, R and p 
are known from experiments. The LHS can be related 
approximately to p/K + and tt + /K + : 



P 

K+ 



N 
K' 



7T+ 

~K+ 



2U 



(21) 



The experimental values of these ratios at SPS are 119 
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= 4.76. 



The values of p is 



p = 0.07. 



(22) 



(23) 



With these experimental inputs there should be no diffi- 
culty in satisfying Eq. (^) and ( |2(i| ) by varying the quark 
numbers. 

However, in our approach the quark numbers must fit 
into our scheme of quark enhancement via gluon con- 
version. Moreover, there is the issue of what precisely is 
the central region where the experimental numbers of the 
hadron ratios are measured. Clearly, the valence to sea 
quark ratio depends on the region of small x considered. 
The experiments do not have a common and unique def- 
inition of the central rapidity region. For our analysis 
in the following we define the central rapidity region to 
correspond to a value xq, which depends on y/s. In fact, 
when we make prediction later on for RHIC energies we 
shall use the relation 



x Q = (s /s) 1/2 . 



(24) 



For now at SPS we use xq as an adjustable parameter. 
The important physics input is that for x < xq, which is 
what Feynman called the "wee" x region, we assume that 
all partons distributions are constant (consistent with the 
notion of saturation) so that the quark number ratios in 
the wee region, whatever the flavor, can be determined 
by computing the ratios of the corresponding quark dis- 
tributions at x = xq. 

In the calculation of the quark distributions after gluon 
conversion, we shall do it in a simplified way for the cen- 
tral region, different from how it has been treated in the 
fragmentation region Jl||] . The reasons are because firstly 
we need not distinguish u and d types quarks in the cen- 
tral region and secondly we need not track the x depen- 
dences. It is important to first refer all quark numbers to 
those given by CTEQ at Q 2 = 1 GeV and x = x - That 
is our baseline, from which we discuss enhancement in 
the following. Now, CTEQ gives distributions, not num- 
ber of partons. For example, u(xq) is the probability of 
having a u quark at x = xq. In the following we use go 
to denote the number of u and d quarks in x < xq, and 
equate it to u(xq) + d{xo), multiplied by a factor that is 
proportional to the relevant phase space volume. Such 
a factor will cancel later upon taking the ratio of quark 
numbers, so it will not appear explicitly in any of the 
expressions for parton numbers below. Similarly, we use 
so and go to denote the number of s quark and gluons in 
the region x < xq, but identified with s(xq) and g(xo), 
respectively, of the CTEQ distributions. 



III. GLUON CONVERSION 

Before gluon conversion we have q v valence quarks, 2q~o 
non-strange sea quarks, and sq + sq strange quarks. In 



the case of hadronic collisions, it has been shown that 
the inclusive distributions of produced hadrons can be 
reproduced without any free parameters, if the gluons are 
completely converted to non-strange sea quarks before 
hadronization through recombination p5[ . Now, in the 
case of AA collisions we must consider the conversion of 
gluons to strange quarks in addition to the non-strange 
quarks because of Pauli blocking in the light sector. We 
use 7 to denote the fraction in the strange sector. That is, 
the number of converted strange and non-strange quarks, 
labeled with subscript c, are 



75o, 



q c = (1-7)30, 



(25) 



with the corresponding antiquarks s c and q c being equal 
in number, respectively. Thus after conversion we have 



q = q v + qo + q c , 

S = Sq + S c . 



(26) 
(27) 



The quark and gluon distributions at xq can be either 
obtained from the graphs posted by CTEQ4LQ [Q, or 
determined numerically from the analytic formulas given 
in Ref. [§2j. We use the latter to fix q v ,qo, sq and go for 
every xq, while q c and s c depend on 7. Hence, we have 
two free parameters, xq and 7, to fit the data through 
the use of Eos. . 

From Eq. ( ]l9| ) one gets q v /s = 3.86. Using that in j20| ) 
yields q/s — 7.53, whereupon one obtains 



l + q v /q 



0.66. 



From 



j) and ( |27|) we have 

s = s + 7 g = q v I 3.86, 
Qo + (1- 7).9o = 7.53 s; 



together they give 



qo + so + go = 2.21 q v . 



(28) 



(29) 
(30) 



(31) 



This is an equation that depends on CTEQ distributions 
only, so we can solve for the value of xq . The result then 
determines also the values of q v , sq and go, which, when 
used in Eq. d2^), fix 7. The process yields 



xo = 0.135, 
7 = 0.08. 



(32) 
(33) 



The value of xo is reasonable, but the value of 7 seems 
surprisingly low, since 8% conversion from the gluons 
seems insufficient to justify the notion of SE. 



IV. STRANGENESS ENHANCEMENT 

To appreciate the value of 7 found above, let us exam- 
ine the quark distributions at xo before gluon conversion. 
Our solution of Eq. (31) gives 



3 



Xq q v — 0.462, x Q q = 0.118, 
xq sq — 0.052, xq go = 0.85. 



(34) 



Thus from Eq. ( p5| ) we have xq s c = 0.068. Comparing s c 
with so 7 wc scc that the strangeness enhancement factor 
E s at the quark level is 



E s = — 

so 



1 



— = 2.3. 

so 



(35) 



This indicates quite an appreciable amount of increase of 
the strange quarks, qualitatively consistent with the hy- 
peron enhancement. The point is that there are so many 
gluons that an 8% conversion significantly enhances the 
strangeness content. The remaining 92% conversion to 
q c should be compared to 100% conversion in the case 
of hadronic collisions fll5|| . Let us call the light quark 
population in the sea after 100% conversion q±, i.e., 



9i = 9o + .9o- 



(36) 



Then the change in the sea from pp to AA collisions can 
be characterized by the ratio B: 



D 



<]_ 



go + (1 - 7) 9o 
9o + .9o 



0.94. 



(37) 



This may be regarded as a numerical factor quantifying 
Pauli blocking in the light quark sector. Note that it 
is less than one by only a small amount, but enough to 
boost E s from one by more than a factor of two. 

The extension of this consideration to RHIC energies 
is straightforward. We first use Eq. (|24|) to determine sq 
from the values of xq at SPS. Setting y/s = 17 GeV, we 
obtain y/s Q = 2.3 GeV. Now, holding sq fixed, we have 
the corresponding xq value (call it x' ) at y/s = 130 GeV 
to be 



x' = 0.0177. 



(38) 



The values of 9„, 90, Sq and go at x' are (from CTEQ) 



x Q q v 



0.149, x' q Q = 0.184, 
x' s = 0.089, x o5o = 1.229. 



(39) 



Note that q' v is much smaller than q v , as expected, so 
that q'0/l'o = 0.55, even before gluon conversion. Thus 
we expect antiparticle/particle ratios to be much closer 
to one. 

As before, we need the experimental inputs at RHIC. 
From Refs. ||^|| we have at y/i = 130 GeV 

r = 1.136, R = 0.77, p = 0.64. (40) 

So we get from Eq. ( [l3| ) 

re = 0.628, re = 0.713. (41) 

Assuming that 7 remains constant, we now can calculate 
the quark ratios 



%=Q. 
s' 



—7 = 7.06, 



(42) 



which, when used in (|19|) and (|20|), enable us to calculate 
the hadron ratios. As a consequence, we obtain 



P 
K+ 



= 0.71, 



K- 



= 0.18. 



(43) 



The latter, compared to the value, 0.21, at SPS, is a 14% 
decrease and agrees well with the data at RHIC [ p6| , 
which shows K + /tt + = 0.176 ± 0.004. For the former 
we find indirect confirmation from the following ratios 
reported by STAR ||0 for Vs = 130 GeV: 



P 
P 



= 0.65 ±0.07, 



P 



7T 



= 0.08 ±0.01, 



— = 0.149 ±0.02, — = 0.88 ±0.05. 

7T K + 

These numbers can be used to imply 



=0.73 ±0.1, 



(44) 



which agrees well with the calculated number in (|43j). 
These results give support to our assumption that 7 is 
constant when the energy is increased and to our proce- 
dure of treating the energy dependence. 

The SE factor becomes at RHIC (y/s = 130 GeV) 



E, 



1 



2.1. 



(45) 



Although the gluon density increases by 45% as Xq de- 
creases to x' Q , the s quark density in the quiescent sea 
increases by even more, so the net enhacement factor 
E s decreases slightly. This small decrease is in agree- 
ment with that of the statistical model j2q], although 
the physics is totally different. The Pauli blocking factor 
becomes 

B = 0.93, (46) 
which is essentially unchanged from Eq. (|37j). 

V. CONCLUSION 

Since we have left out the multi-strange hyperons from 
our consideration, we cannot expect the numbers cal- 
culated to be accurate. Moreover, the necessity to use 
such experimental inputs as r, R, p to determine k and 7 
renders the approach highly phenomenological, far from 
first principles. However, the basic attributes of this line 
of study are to use the parton model (and the distribu- 
tions of CTEQ) as the basis for the investigation of par- 
ticle ratios in nuclear collisions at the quark level, and to 
use simple linear relations, Eqs. based on quark 

counting as the only constraints among the strange and 
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non-strange hadrons. We have found consistency within 
this simple approach, and can successfully describe the 
energy dependence. We have not assumed thermal equi- 
librium, nor relied on the statistical model. We have 
also deliberately avoided treating mesons and baryons as 
products of quark densities, as have been attempted in 
Refs. [|]](|, since they lead to either s ^ s or undeter- 
mined constants. 

As we have stated at the outset, it is not our aim to 
predict particle ratios. We have used the experimental 
values of the ratios to lead us to the determination of the 
SE factor, E s , and the Pauli blocking factor, B, defined at 
the quark level. In so doing we have gained some insight 
into how the enhancement mechanism works through the 
process of gluon conversion. We have further learned that 
a slight suppression of the conversion into the non-strange 
sector gives rise to a substantial increase in the strange 
sector. Such a small change from hadronic to nuclear 
collisions makes strangeness enhancement an unreliable 
signature for the formation of quark-gluon plasma. 
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